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Notation 


Symbol 

A  Aspect  ratio,  =  L/D 

Cp  Constant  pressure  specific  heat 

D  Width  of  enclosure 

g  Acceleration  due  to  gravity 

g  Dimensional  constant  in  Newton's  second  law 
o 

Grashof  number,  =  gg( 8jj-0^,)D3 /v2 
h  Local  heat  transfer  coefficient,  =  q/e^-e^ 
k  Thermal  conductivity 

L  Height  of  enclosure 

Nu  Local  Nusselt  number,  =  g— 

Nu  Mean  value  of  Nusselt  number 

p  Local  dynamic  pressure 

__  T~\  2. 

P  Dimensionless  pressure,  = 

Pr  Prandtl  number,  =  gcviCp/k 

q  Heat  flux  at  hot  wall 

S  Inclination  angle 

t  Time 

u  Velocity  in  x-direction 

U  Dimensionless  velocity  in  X-direction,  =  uD 
v  Velocity  in  y-direction 

V  Dimensionless  velocity  in  Y-direction,  =  vd 
x  Distance  along  hot  wall 

X  Dimensionless  distance  along  hot  wall,  =  x/ 


Symbol 


y  Distance  from  hot  wall 

Y  Dimensionless  distance  from  hot  wall,  =  y/D 
Greek  Letters 


At 

C 

0 

0 

v 

v 

p 

T 

4» 


pc 


a  Molecular  thermal  diffusivity,  = 

0  Volume  coefficient  of  thermal  expansion 
AX  Grid  spacing  in  X  direction 
AY  Grid  spacing  in  Y  direction 
Time  increment 


Dimensionless  vorticity,  =  +  -gy) 

Temperature  (0  ,  and  0p  refer  to  the  temperatures  at 
the  hot  and  cold  walls  respectively ) . 

Dimensionless  temperature,  =  ( 0-0..)/ ( 

11  U  n 

Kinematic  viscosity 

Viscosity 

Density 

Dimensionless  time,  =  tv/D2 

Dimensionless  stream  function,  such  that  U=3'}'/3Y  and 

v=  -at/ax 


a)  Relaxation  parameter 
<J>  Slat  angle 
Subscripts 

i,j  Space  grid  point  indices  in  X  and  Y  directions 

opt  Value  of  relaxation  parameter  giving  fastest  conver¬ 
gence. 

w  Wall  grid  point 

Superscripts 

n  Time  index 
m  Iternation  number 
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Abstract 


A  numerical  investigation  was  conducted  on  two-dimen¬ 
sional  natural  convection  within  inclined  rectangular  en¬ 
closures  partitioned  into  45  degree  triangular  cells.  The 
time  dependent  governing  equations,  vorticity,  energy,  and 
stream  function,  were  solved  by  an  ADI  method  and  a  Gauss- 
Seidel  SOR  technique.  The  numerical  procedure  was  validat¬ 
ed  for  rectangular  enclosures,  then  modified  for  triangular 
cells.  Heat  transfer  coefficients  were  determined  for  an 
inclined  square  enclosure  with  a  diagonal  partition  for 
Grashof  numbers  less  than  2.xl05  and  inclination  angles 
between  10  degrees  and  90  degrees.  These  results  show  a 
diagonal  partition  reduces  the  heat  transferred  by  natural 
convection  across  an  inclined  square  enclosure  by  more  than 
50%. 
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NUMERICAL  SOLUTION  OF  NATURAL  CONVECTION 
IN  AN  INCLINED  RECTANGULAR  CAVITY 
WITH  PARTITIONS 

I .  Introduction 

Background 

A  major  factor  in  the  design  of  solar  collectors  is 
the  reduction  of  natural  convection  heat  losses  within  the 
collector  enclosure.  Traditional  designs  have  suppressed 
natural  convection  by  developing  large  aspect  ratio  enclo¬ 
sures  or  designing  a  small  aspect  ratio  honeycomb  structure 
to  separate  the  collector  plates  (Ref.l).  However,  recent 
experimental  studies  by  Meyer,  el.al.  (Ref. 2),  and  Holland 
(Ref. 3)  have  demonstrated  the  effectiveness  of  inclined 
slats  (partitions)  in  reducing  natural  convection  heat  loss¬ 
es.  The  slats  in  these  studies  were  positioned  to  form  a 
series  of  moderate  aspect  ratio  parallelogram  enclosures. 
Meyer,  et.al.,  also  concluded  that  slats  oriented  downward 
from  the  hot  plate  resulted  in  a  reduction  of  convective 
heat  transfer  with  the  minimum  heat  loss  occurring  at  a 
slat  angle  of  45  degrees. 

Objective 

In  the  present  study,  the  influence  of  slats  in 


reducing  natural  convection  will  be  determined  for  non¬ 
parallel  slat  arrangements.  Specifically,  the  slats  will 
be  oriented  so  as  to  partition  a  moderate  aspect  ratio  rec¬ 
tangular  enclosure  into  a  series  of  triangular  regions  (see 
Figure  1).  The  temperature  distribution  and  the  cell  heat 
transfer  coefficients  will  then  be  determined  theoretical¬ 
ly  by  a  numerical  procedure  which  will  also  be  developed 
for  this  investigation.  The  heat  transfer  coefficient  for 
the  triangular  cells  will  then  be  compared  to  results  ob¬ 
tained  for  similar  rectangular  enclosures  with  and  without 
slats.  From  these  results  the  effectiveness  of  triangular 
slat  arrangements  will  be  determined. 

The  first  step  in  this  investigation,  however,  will  be 
to  adapt  a  numerical  procedure  which  is  capable  of  solving 
the  nonlinear  partial  differential  equations  which  describe 
natural  convection  in  an  inclined  enclosure.  The  implicit 
finite  difference  computational  scheme  developed  by  Wilkes 
and  Churchill  (Ref. 4)  is  just  such  a  procedure.  Since  the 
Wilkes  and  Churchill  numerical  procedure  was  developed  for 
a  rectangular  cavity  with  one  vertical  wall  hot  and  the 
other  cold,  the  method  will  be  first  modified  to  describe 
an  inclined  rectangular  cavity,  then  modified  again  to  ac¬ 
count  for  the  addition  of  partitions  in  the  rectangular  en¬ 
closure. 
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Figure  1.  Triangular  Partitioned  Enclosure. 


Scope 


In  this  study,  natural  convection  heat  transfer  coeffi¬ 
cients  were  evaluated  for  moderate  aspect  ratio  (A=1.0  to  3.0) 
rectangular  enclosures  with  thin  partitions.  The  partitions 
were  oriented  at  ±45  degress  to  the  hot  wall.  Thermal  con¬ 
duction  was  assumed  to  be  infinite  across  the  partitions  and 
zero  along  the  partitions.  Additionally,  the  steady  state 
temperature  distribution  along  a  partition  was  considered 
linear.  Finally,  since  only  heat  transferred  due  to  nat¬ 
ural  convection  was  to  be  evaluated,  thermal  radiation  heat 
transfer  within  the  enclosure  was  neglected. 

To  follow  the  development  of  the  present  investigation, 
this  report  is  organized  into  five  additional  chapters. 

The  time  dependent  governing  equations  which  describe  nat¬ 
ural  convection  in  inclined  rectangular  cavities  are  devel¬ 
oped  in  Chapter  II.  In  Chapter  III,  a  numerical  procedure 
is  developed  to  evaluate  the  governing  equations.  A  dis¬ 
cussion  of  the  numerical  results  is  presented  in  Chapter  IV. 
Finally,  the  significant  conclusions  are  restated,  and  recom¬ 
mendations  are  made  for  further  studies  in  Chapter  V  and  VI, 
respectively. 
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II.  Mathematical  Model 


In  this  section  the  governing  equations  describing  nat¬ 
ural  convection  in  a  2-D  inclined  rectangular  cavity  will  be 
developed.  The  equations  and  boundary  conditions  will  then 
be  simplified  and  rewritten  in  nondimensionalized  form.  Fin¬ 
ally,  an  expression  for  local  and  average  heat  transfer  co¬ 
efficient  will  be  derived. 

Rectangular  Cavity 

Consider  the  inclined,  rectangular  cavity  in  Figure  2. 
The  left  and  right  walls  are  at  constant  temperatures  e„,  and 

n 

0 ,  respectively,  and  the  other  side  walls  are  either  insul- 
ated  or  have  a  linear  temperature  distribution.  Initially, 
the  fluid  in  the  cavity  is  motionless,  and  at  a  constant 
temperature  equal  to  the  average  temperature  of  the  hot  and 
cold  walls.  The  fluid  thermodynamic  and  transport  proper¬ 
ties  except  density  are  considered  constant  and  independent 
of  temperature.  Under  the  Boussinesq  Approximation,  density 
changes  are  expressed  only  as  a  function  of  temperature  dif¬ 
ferences,  and  density  is  assumed  constant  except  in  the  buo¬ 
yant  force  terms. 

Under  these  conditions,  the  fluid  in  an  inclined  rec¬ 
tangular  cavity  is  described  by  the  following  equations: 
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Conservation  of  Mass 


3u  3v  _ 
3x  3y 


0 


(1) 


Conservation  of  Momentum 


p  3u  P  u3u  v3u  _ 

gc  3t  gc  ax  gc  ay 


_  l£  _  JB£ 


ax  gc 


sin 


p  3v  (  P  u3v  +  jP_  v3v  _ 
gc  3 1  gc  3x  gc  ay 


ap 

ay 


_  JS 

gc 


(3) 


Conservation  of  Energy 


30 

3t 


+ 


a( 


32e 

3X2" 


32  0\ 

ayv 


(4) 


Coefficient  of  Volumetric  Expansion 


0 


1  3p 


P  o-P 

p(0-eo) 


(5) 


The  associated  boundary  and  initial  conditions  are  as 
follows : 


0  (  x ,  0 ,  t ) 

(6) 

0 (x, D, t )  =  ec 

(7) 
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Linear  Temperature  Side  Walls 


9(0, y,t)  =  0H+(0H-0c)y/D 


0(L, y , t )  =  0H+(0H-0c)y/D 


Adiabatic  Side  Walls 


(0,y,t)  = 


||(L,y,t)  -  0 


Cavity  Wall  No-Slip  Condition 


u(x,0,t)  =  u(x,D,t)  =  0 
v(0,y,t)  =  v(L,y,t)  =  0 


Initial  Conditions 


u(x,y , 0)  = 


v(x,y,0)  = 
0(x,y,O)  = 


0H+0C 


(ID 


Since  the  pressure  gradient  in  the  cavity  results  in  part 
from  the  change  in  elevation  of  the  fluid,  it  is  conven¬ 
ient  to  derive  an  expression  for  the  total  pressure  grad¬ 
ient  in  terms  of  the  dynamic  pressure  gradient  and  the 
change  in  pressure  due  to  the  weight  of  the  fluid. 
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p  =  p0_  |£±^(Sin  S)x+(cos  S)y] 
sc 


(12) 


Differentiating  this  expression  yields 

If  =  -  f-  p°  Sin<s>  <13) 

o 

If  =  ~fy  “  f“  Po  cos(S>  (I4) 

o 


Substituting  from  Equations  (13), (14)  and  (  5)  into  equations 
(2)  and  (3)  and  simplifying 


3u 

3t 


+ 


£c  apo.  _ 
p  3x 


8(9o-e 


)(15) 


3v 

3t 


+ 


£c  3p  0 
p  3y 


8(80-9)g  cos 


+ 


)(16) 


Equations  (1),  (15),  (16)  and  (4),  and  the  associated 


boundary  and  initial  conditions  can  be  restated  in  nondimen- 


sional  form.  The  nondimensional  parameters  are  defined  in 
the  List  of  Symbols. 


3U  3V 
3X  3Y 


0 


(17) 
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■x 


40  * 


M  +  ui2  +  VM  = 

3r  U3X  V3Y 


_  3P  _  G  0  sirWSU92u  +  9'U 

3X  k  0  sm(b)+_r  +  —y 


IX  +  u9V  V3V  = 

3t  U3X  3Y 


3P  _  _  .  / o \ . 3 2  V  3  2  V 

3Y  "  GR  0  sln(S)  3X2-  +  ay2- 


(19) 


30  n30  30  1  ra2e  3 2  0 

3t  3X  v  3Y  P  3Xr  W* 

r 


(20) 


0(X, 0, r )  =  0 


0(X,1,T)  =  1.0 


0(0, Y,T)  =  0(A,Y,t)  =  Y 


H(0,Y,t)  =  |§(A,Y,t)  =  0 


(21) 


U(X, 0, t )  =  U(X,1,T)=  V( A , Y , t )  =  V( 0, Y, x )  =  0 


U(X, Y, 0)  =  V( X , Y , 0 )  =  0 


0(X,Y,O)  =  0.5 


The  pressure  terms  in  Equations  (18)  and  (19)  can  now 
be  eliminated  by  differentiating  Equation  (18)  and  (19)  with 
respect  to  Y  and  X  ,  respectively,  subtracting  and  simpli¬ 
fying  the  resulting  expression  with  Equation  (17).  This  pro¬ 
cess  yields  Equation  (22). 


_lf3U  3Vl+II[3U  3Vl,vf3U  3V] 
3t|_3Y  "  3XJ  UL3Y  ~  3Xj  [' 3Y  "  3xJ 

-  Gr[cos(S>H  -  3in<s,|f]+  £[%  -  ||] 


+  ^ 

+  TP 


f3U  3Vl 
[_3 Y  “  3XJ 


(22) 


If  a  vorticity  function  (?)  is  now  defined  as 


3U  3V 
~Q  ~  3Y  3X 


(23) 


and  the  stream  function  (M')  is  introduced,  where 


(24) 


(25) 


Equation  (22)  can  be  rewritten  in  terms  of  vorticity,  and 
Equation  (23)  in  terms  of  the  stream  function. 


ii  +uh  +vi£ 

3t  3X  3Y 


;[sin(S)||  - 


cos(S)-! 


.  3 2  C  .  32C 

Tx2”  TP 


_  _  3Y,  3? 
~  7X  3Y 


(26) 


The  governing  equations  and  boundary  conditions  which  de¬ 
scribe  natural  convection  in  an  inclined  rectangular  cavity 
can  now  written  in  their  final  form. 

Vorticity  Equation 

3t  U3X  V3Y  GR[Sin(S)3Y  “  Cos(S)3X_ 


32?  3  2  £ 

+  — Z-  +  — 4- 
3X7  3Y7 

(28) 

Energy  Equation 

il&  +-U—  +V—  = 
3t  U  3X  3Y 

1  [3Z0  320l 

PR  [3X7  W7] 

(29) 

Stream  Function  Equation 

-c  = 

3  2  'F  3  2  Y 

3X7  3Y7 

(30) 

Velocity  Component  Equations 

u  = 

U  3Y 

(31) 

V  —31 

3X 

(32) 
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Boundary  Conditions 


U(0,Y,t)  =  U( A , Y , t )  =  0 
V(X,0,x)  =  V(X,l,x)  =  0 

0,  Y,  x  )  =  Mf(A,Y,T)  =  Y(X,  0,  t  )  =  Y(X,  1,  t  )  •  0 
0(X,O,x)  «  0 
0(X,1,T)  -  1.0 
Adiabatic  Side  Walls 

||(0.Y,t)  =  ||(A,Y,t)  =  0 

Linear  Temperature  Side  Walls 

0(0, Y,T)  *  0( A , Y , t )  =  Y 

Initial  Conditions: 

U(X, Y , 0 )  =  V(X, Y, 0)  =  0 

nX,Y,0)  =  C(X,Y,0)  =  0 

0( X , Y , 0 )  =  0.5 

Heat  Transfer  Coeffi c ien t_s 

Local  and  average  heat  transfer  coefficients  are  ex' 
pressed  as  non-dimensional  Nusselt  ..umbers. 


Nu  -  hD/k 


Local  Nusselt  Number 


The  convection  heat  transfer  coefficient  (h)  at  the 
hot  wall  is  defined  as  follows: 


h 


-koe/3y)waU 


0H  '  9C 


(38) 


Introducing  the  definition  of  Nu ,  Equation  (37),  0  ,  and 

Y  ,  Equation  (38)  becomes 


Nu 


(ifi) 

v  9 Y ''wall 


(39) 


Average  Nusselt  Number 

To  obtain  an  average  Nusselt  Number,  Nu,  along  the  hot 
wall,  the  local  Nusselt  Numbers,  Nu,  are  integrated  over  the 
length  of  the  hot  wall. 


(40) 


Partitioned  Cavity 

A  review  of  the  mathematical  model  of  natural  convec¬ 
tion  in  an  inclined  rectangular  enclosure  shows  the  govern¬ 
ing  Equations  (28),  (29),  (30),  (31)  and  (32)  are  indepen¬ 
dent  of  the  enclosure  geometry.  Therefore,  these  equations 
can  be  applied  to  triangular  enclosures,  and  only  the  boun¬ 
dary  conditions  need  to  be  changed. 
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The  partitions  will  be  placed  in  the  enclosure  at 
<J>=±45  degrees  with  respect  to  the  hot  wall,  and  also  all 
partition  steady  state  temperature  distributions  will  be 
linear.  Therefore,  applying  these  assumptions,  the  boun¬ 
dary  conditions  for  the  partitions  can  now  be  written. 

<J>  =  45° 

0( 1. 0-Y, Y, x )  =  Y  (41) 

Y(1.0-Y,Y,t)  =  0 


U(1.0-Y,Y,T)sin<}>+V(1.0-Y,Y,T)sin<J>  =  0 

4>  =  -45 
0(Y,Y,t)  =  Y 
Y(Y,Y,t)  =  0 
U(Y,Y,T)sin<J)+V(Y,Y,T)sin<j>  =  0 


(42) 


Now  the  partition  boundary  conditions  can  be  used  in  com¬ 
bination  with  those  for  the  rectangular  cavity  to  complete¬ 
ly  describe  the  boundaries  of  an  inclined  rectangular  en¬ 
closure  with  partitions. 
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III.  Numerical  Procedure 


Rectangular  Enclosure 

This  section  contains  a  description  of  the  procedure 
that  was  used  to  solve  the  system  of  equations  developed 
in  the  previous  section.  This  procedure  was  adapted  from 
a  finite  difference  computational  procedure  developed  by 
Wilkes  and  Churchill  (Ref. 4)  to  study  natural  convection 
in  a  rectangular  enclosure  with  one  vertical  wall  heated 
and  the  other  cooled.  This  procedure  is  illustrated  by  a 
flow  chart  in  Figure  3.  Because  the  procedures  are  so  sim¬ 
ilar,  tests  of  stability  and  convergence  were  not  rigorous¬ 
ly  carried  out  in  this  study. 

The  geometry  of  the  rectangular  enclosure  and  the  fin¬ 
ite  difference  nomenclature  are  shown  in  Figure  4.  The 
mesh  spacing  in  the  X-direction  is  AX  and  the  Y-direction 
is  AY  .  Subscripts  (i,j)  are  associated  with  each  mesh 
point,  so  that  the  space  variables  may  be  expressed  as 
X .  =( i-1 )AX  ,  for  i=l,2,...,I  and  Y.=(j-1)AY  ,  for  j= 

1,2,3,...,J.  Time  is  segmented  into  equal  intervals  At 
so  that  the  nondimensional  time  x  is  x=nAx  for  n=0,l,2. 

The  notation  4)1?  .  denotes  the  values  of  the  variable 
^  >  J 

at  mesh  point  (i,j)  at  time  level  n 

The  first  equation  to  be  solved  from  nAx  to  (n+l)AT 
is  the  energy  equation  (29),  which  is  parabolic  in  time, 
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Figure  3.  Numerical  Procedure  Sequence  of 
Solutions  of  the  Governing  Equations. 


17 


_  r* •» 


- 1  ‘lirrnBi  t> i  * 


Figure  3.  Continuation  of  Numerical  Procedure  Sequence 
of  Solution  of  the  Governing  Equations. 

nonlinear  and  of  second  order.  The  solution  to  this  equa¬ 
tion  is  advanced  one  time  level  by  using  an  implicit  alter¬ 
nating  direction  (ADI)  technique  developed  by  Peaceman- 
Rachford  (Ref.7). 

To  apply  the  ADI  method  to  the  energy  equation  (29)  re¬ 
quires  the  coefficient  velocities  U  and  V  to  be  held  con¬ 
stant  at  any  grid  point  over  a  time  step.  Now  the  partial 
derivatives  in  Equation  (29)  can  be  approximated  by  finite- 
differences.  (Note:  All  space  derivatives  are  central  dif¬ 


ferences  .  ) 


T=(n+£)AT 


Equations  (43)  and  (44)  are  implicit  in  X  and  Y  direc¬ 
tions,  respectively,  and  when  applied  to  every  point  in  a 
row  or  column,  as  the  case  may  be,  yield  a  tridiagonal  matrix 
in  the  unknown  temperatures  or  .  In  this  study, 

a  tridiagonal  algorithm  adapted  by  Roache  (Ref. 7)  was  used  to 


solve  the  tridiagonal  matrices. 

The  vorticity  equation  (28}  which  has  characteristics 
similar  to  the  energy  equation^  can  now  be  solved  by  the  ADI 
method  for  new  values  of  vorticity  at  the  n  inner  grid 
points.  The  vorticity  at  the  boundary  points  will  remain 
at  their  old  nAx  values  during  this  calculation,  and  will 
be  advanced  later.  Because  there  is  a  temperature  gradient 
in  the  vorticity  equation  (28),  the  new  (n+l)Ax  tempera¬ 
ture  field  will  be  used  to  evaluate  the  temperature  gradient 
at  each  grid  point.  The  temperature  gradient  will  then  be 
held  constant  throughout  the  calculation.  The  finite  differ¬ 
ence  representations  of  the  vorticity  are 

T~(n+£)AT 
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x=(n+l)Ar 


n+1  n+i 

Zhl—hA 

Ax/2 


+  un  . 

i»  J 


~n+i  n+£ 

’n+1  n+1 

^i+l,j  "  H-l ,j 

+vn 

H,j  +  1  "  H.j-l 

2AX 

'  V  .  . 

1  >  J 

2AY 

=  Gr 


sin(S) 


r0n  +  1  - 

2AY 


-  cos(S)l 


’0n+^  -  0n+^ 
l+l ,  j 

2AX 


rn+h 

6i+l.d' 

-2t,n+i.  + 
i,  J 

n+i 

+ 

\n+1  -2rn+1  +  rn+1 

?i>j+l  ^i,j  ^ i , j - 1 

AX2 

AY2 

(46) 


Equations  (45)  and  (46)  are  also  implicit  in  X  and  Y 
directions,  respectively,  and  when  applied  to  every  point  in 
a  row  or  column,  as  the  case  may  be,  also  yield  a  tridiagonal 
matrix.  These  matrices  are  solved  by  the  same  algorithm  that 
was  used  for  solving  the  energy  equation. 

With  the  new  interior  vorticity  field  calculated,  the 
method  of  successive  over-relaxation  is  used  to  solve  the 
stream  function  equation  for  the  new  stream  function  field. 
The  expression  for  the  m+1  iteration  of  the  stream  func¬ 
tion  at  a  point  is  as  follows: 
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TABLE 

I 

OPTIMUM  VALUES  OF  RELAXATION  PARAMETER 

FOR  USE  IN  EQ . (47 ) 
(Ref .4) 

No.  of  grid 

spacings : 

X  Direction 

Y  Direction 

> 

X 

II 

> 

w0PT 

10 

10 

0.2 

1.58 

20 

20 

0.05 

1.83 

20 

10 

0.1 

1.70 

30 

10 

0.1 

1.73 

Wilkes  and  Churchill  (Ref. 4)  determined  optimum  values 
of  relaxation  parameters  for  some  representative  grid  sizes 
as  seen  in  Table  I.  Using  these  values,  Wilkes  and  Church¬ 
ill  experienced  good  convergence  at  each  grid  point  with 
about  twenty  five  iterations.  In  this  study  twenty  five 
iterations  were  performed  for  each  grid  point  using  the  op¬ 
timum  value  of  the  relaxation  parameters  in  Table  I. 

The  new  wall  vorticities  can  now  be  determined  from  the 
new  stream  function  field.  In  this  study,  two  wall  vorticity 
expressions  were  used  depending  upon  the  Grashof  Number  of 
the  fluid.  Wilkes  and  Churchill  results  indicated  the  first 
expression  will  cause  instabilities  above  GR=200,000. 

n+1  -  V?+l 

5  = - =*= - “  when  G  <200,000  (48) 

i , 0  2( AY)2 
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when  G^>200,000 


(49) 


n+1 

2)  C 

i,0 


24'n+4 

-  1  >  y 

(AY V 


Both  expressions  were  derived  by  Wilkes  and  Churchill  (Ref. 4) 
The  new  velocity  fields  U  and  V  can  now  be  calcul¬ 
ated  from  the  new  stream  function  field.  In  this  study  sec¬ 
ond  order  space  central  finite  difference  approximations  were 
used  for  the  new  velocities. 


H* .  .  -2'}' .  .+¥.  .  , 

nn+1  -  1..1  +  1  i.J  1..1-1 

i,j  ~  2AY 


(50) 


ur  —2^  +4' 

vn+l  =  -  i+1.J  i  .3  i-1.,1  (5 

1,J  2AX 

These  expressions  differ  from  those  of  Wilkes  and  Churchill 
who  used  a  fourth  order  finite  difference  approximation. 

Now  using  the  newest  temperature  field,  the  local  and 
average  Nusselt  Numbers  can  also  be  calculated.  The  deriv¬ 
ative  in  Equation  (39)  is  expressed  by  a  second  order  for¬ 
ward  finite  difference  relation. 


Nu 


ei,3-40i,2+30i,l 

2AY 


(52) 


The  average  Nusselt  Number  is  obtained  by  integration  using 
the  Trapezoidal  rule. 
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Partitioned  Cavity 

In  the  previous  chapter,  the  mathematical  model  of  nat¬ 
ural  convection  in  an  inclined  rectangular  enclosure  with 
partitions  was  developed.  The  governing  equations  were  de¬ 
termined  to  be  unchanged.  However,  new  boundary  conditions 
were  necessary  to  describe  the  partitioned  enclosure.  In 
this  section,  the  numerical  procedure  previously  developed 
will  be  modified  to  account  for  45  degree  partitions.  And 
finally,  a  method  of  solution  of  a  partitioned  rectangular 
enclosure  will  be  discussed. 

Numerical  Procedure  Modification.  The  two  significant 
modifications  made  to  the  inclined  rectangular  enclosure 
numerical  procedure  were;  1)  a  new  wall  vorticity  model  for 
the  45  degree  partition,  and  2)  local  heat  transfer  coeffi¬ 
cients  for  the  partition.  While  the  wall  vorticity  models 
used  in  the  rectangular  enclosure  are  accurate  for  grid 
points  along  the  mesh,  the  following  expression  (see  Ref. 7 
page  144)  provides  a  more  accurate  evaluation  of  the  wall 
vorticity  of  a  sloping  wall  (<{>=45  degrees)  with  AX=AY. 


^iw, jw 


2(4'.  ,  .  +4».  .  -24' .  .  ) 

v  1W-1,JW  iw,jw+l  1W, Jff' 

AX 


(53) 


This  expression  was  used  to  evaluate  the  wall  vorticity  of 
all  partitions. 

Finally,  the  local  heat  transfer  equation  was  modified 


to  determine  the  heat  transfer  coefficient  across  a 


partition.  First  and  second  order  finite  difference  schemes 


were  sufficient  for  all  grid  points  along  the  partition  ex¬ 
cept  the  end  points.  At  the  end  points,  the  heat  transfer 
coefficient  across  the  partition  was  determined  to  be  equal 
to  the  heat  transferred  along  the  partition.  Details  of  the 
heat  transfer  coefficients  expressions  can  be  found  in  Ap¬ 
pendix  B. 

Solution  Method.  Since  the  natural  convection  govern¬ 
ing  equations  did  not  change  with  the  addition  of  parti¬ 
tions,  each  triangular  cell  can  be  independently  solved  by 
assuming  common  boundary  conditions  for  the  neighboring  tri¬ 
angles.  However,  it  is  then  necessary  to  compute  an  energy 
balance  at  each  partition  grid  point.  If  the  partition  boun¬ 
dary  conditions  do  not  balance,  it  is  necessary  to  iterate 
the  triangular  cell  partition  boundary  conditions  until  the 
heat  transfer  coefficients  and  temperatures  match.  This  it¬ 
erative  procedure,  however,  was  simplified  in  this  study  by 
the  thermal  characteristic  of  the  partition  model. 

In  this  study,  the  partition  model  steady  state  temp¬ 
erature  distribution  was  linear,  and  the  partition  was  a 
pure  conductor  across  the  partition.  Because  of  these  ther¬ 
mal  characteristics,  the  temperature  and  heat  transfer  co¬ 
efficients  at  each  partition  grid  point  always  match. 
Therefore,  it  was  decided  to  build  a  composite  model  of  the 
rectangular  enclosure  with  partitions. 


The  composite  model  was  composed  of  individual  tri¬ 
angular  cells  which  were  each  solved  independently.  Then 
the  solutions  were  added  together  to  obtain  the  final  nat¬ 
ural  convection  heat  transfer  results.  A  complete  compu¬ 
ter  program  listing  is  in  Appendix  B. 


IV.  Discussion  of  Results 


The  discussion  of  results  centers  on  two  areas:  deter¬ 
mining  the  validity  and  accuracy  of  the  numerical  procedure 
and  computer  code;  and  determining  the  effect  of  partition¬ 
ing  an  inclined  rectangular  enclosure  into  triangular  re¬ 
gions  on  natural  convection  heat  transfer. 

The  computations  were  conducted  (using  a  CDC  6600  com¬ 
puter)  for  air  cavities  (Pr=.7333)  with  aspect  ratios  from 
1.0  to  3.0  at  heat  fluxes  yielding  values  up  to  2xl05. 

Rectangular  Enclosure  Numerical  Results 

To  determine  the  validity  and  accuracy  of  the  numeri¬ 
cal  procedure,  transient  and  steady  state  heat  transfer  re¬ 
sults  were  obtained  for  rectangular  enclosures  at  various 
inclination  angles,  aspect  ratios,  side  wall  boundary  con¬ 
ditions  and  Grashof  numbers.  These  results  were  then  com¬ 
pared  with  previous  numerical  results  obtained  by  Wilkes 
and  Churchill  (Ref. 4),  Koutsoheras  (Ref. 5),  and  Ozoe ,  et.al. 
(Ref. 6).  A  complete  summary  of  the  computer  runs  is  pre¬ 
sented  in  Appendix  A. 

Vertical  Rectangular  Enclosure.  The  accuracy  of  the 
heat  transfer  solution  was  determined  by  comparing  the  steady- 
state  numerical  solution  of  the  rectangular  enclosure  in  the 
vertical  position  (S=90°) 
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TABLE 

II 

Comparison  of 

Steady  State  Mean  Nusselt  Numbers 

Grashof 

Number 

(Gr) 

Aspect 
Rat  io 
(A) 

Side  wall 

Boundary 

Condition 

Mean  Nusselt 
Number  (Nu) 

Wilkes  & 
Churchill 
Mean  Nus¬ 
selt  no. 
Nu 

Per¬ 
cent 
Diff . 

20000 

2.0 

Insulated 

3.036 

2.874 

5.60 

20000 

2.0 

Insulated 

3.141 

2.992 

4.98 

20000 

3.0 

Insulated 

2.954 

2.825 

4.56 

60000 

1.0 

Insulated 

4.905 

4.793 

2.34 

100000 

1.0 

Insulated 

6.136 

5.512 

11.30 

20000 

1.0 

Linear 

2.061 

2.068 

0.34 

with  heat  transfer  results  attained  by  Wilkes  and  Churchill 
(Ref.4).  A  comparison  of  the  two  results  is  shown  in  Table 

II. 

The  results  in  Table  II  show  that  in  each  comparison, 
the  numerical  procedure  developed  in  this  study  overestimates 
by  five  to  12  percent  the  mean  Nusselt  numbers  achieved  by 
Wilkes  and  Churchill.  This  difference  in  results  could  be 
caused  by  one  or  both  of  two  differences  in  the  numerical 
procedures . 

In  this  investigation,  the  velocity  field  was  calculat¬ 
ed  from  the  new  stream  functions  with  a  second  order  finite 
difference  representation.  On  the  other  hand,  Wilkes  and 
Churchill  used  a  fourth  order  finite  difference  model  to 
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calculate  the  velocity  field.  The  higher  order  model  yields 
more  accurate  velocities  per  time  step.  But,  since  this  in¬ 
vestigation  was  primarily  interested  in  the  steady  state  sol 
ution,  the  improved  accuracy  of  the  transient  solution  vel¬ 
ocity  field  was  not  necessary.  So  it  was  decided  to  use  the 
less  complicated  second  order  finite  difference  representa¬ 
tion. 

A  second  difference  between  the  numerical  methods 
could  be  the  finite  difference  representation  used  to  deter¬ 
mine  the  local  Nusselt  numbers  at  the  hot  wall.  The  wall 
heat  flux  approximation  used  by  Wilkes  and  Churchill  (Ref. 4) 
was  not  specified  in  their  report.  So,  it  was  decided  to 
calculate  the  local  heat  flux  with  a  se  and  order,  forward 
difference  approximation  in  this  study.  The  second  order  ap 
proximation  will  provide  sufficient  accuracy  for  comparison 
and  trend  information,  but  may  represent  the  difference  ob¬ 
served  in  the  heat  transfer  results. 

One  further  comparison  of  the  two  numerical  procedures 
was  to  compare  the  transient  heat  transfer  results.  Figure 
5  is  a  plot  of  this  study's  mean  Nusselt  numbers  vs.  non- 
dimensional  time,  starting  from  the  initial  condition  to 
steady  state.  The  shape  and  indicated  trends  of  the  trans¬ 
ient  solution  curve  is  identical  to  the  plot  achieved  by 
Wilkes  and  Churchill  (Ref. 4).  A  conclusion  that  can  be  de¬ 
duced  from  this  result  is  that  eventhough  the  two  numerical 
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Figure  5.  Transient  Solution  for  the  Vertical  Square  Enclosure 


procedures  produce  slightly  different  values  for  mean  Nus- 
selt  numbers,  the  difference  is  constant  throughout  the  sol¬ 
ution.  Therefore,  the  numerical  procedure  developed  in  this 
study  is  an  adequate  numerical  model  of  transient,  as  well 
as  steady-state,  natural  convection  in  a  vertical  rectan¬ 
gular  enclosure. 

Inclined  Rectangular  Enclosure.  The  previous  discus¬ 
sion  was  limited  to  determining  the  adequacy  of  this  study’s 
numerical  procedure  for  the  vertical  rectangular  enclosure. 
Now  it  is  necessary  to  analyze  the  results  of  the  numerical 
procedure  and  computer  code  over  the  total  planned  operat¬ 
ing  range  of  each  of  the  variables.  Figures  (6),  (7),  and 
(8)  illustrate  the  effects  of  inclination  angle,  aspect  ra¬ 
tio,  Grashof  number,  and  side  wall  boundary  conditions  on 
heat  transfer.  The  following  are  conclusions  which  can  be 
drawn  from  the  figures. 

I.  Effect  of  Side  Wall  Boundary  Conditions.  Figures 
(6)  and  (7)  show  that  the  heat  transfer  at  the  hot  wall  is 
a  strong  function  of  the  side  wall  boundary  conditions.  The 
heat  transferred  at  the  hot  wall  decreases  approximately  35% 
when  the  side  walls  are  changed  from  an  insulated  boundary 
to  a  pure  conductor  across  the  boundary  (linear  temperature 
distribution  along  the  side  wall).  The  effect  of  the  side 
wall  boundary  condition  also  appears  to  be  independent  of 
Grashof  number  and  inclination  angle.  The  results  of  this 
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Grashof  Number,  Gr 

Figure  6.  Effect  of  Side  Wall  Boundary 
Conditions  on  Heat  Transfer. 

effect  are  in  excellent  agreement  with  those  reported  by 
Koutsoheras  (Ref. 5). 

II.  Effect  of  Inclination  Angle,  The  effect  of  inclina¬ 
tion  angle  is  best  illustrated  by  Figure  7.  The  highest  heat 
transfer  takes  place  when  the  enclosure  is  inclined  at  60° 
from  the  horizontal.  However,  the  total  change  in  heat  trans¬ 
fer  due  to  a  change  in  inclination  angle  is  insignificant  when 
compared  to  side  wall  boundary  conditions  and  aspect  ratio 
changes.  This  result  was  previously  observed  by  Koutsoheras 
(Ref. 5)  and  Ozoe  (Ref. 6). 

III.  Effect  of  Aspect  Ratio.  The  effect  of  aspect  ratio 


on  heat  transfer  for  the  case  of  moderate  aspect  ratio  enclo¬ 
sures  with  insulated  side  walls  is  illustrated  in  Figure  8. 


This  figure  shows  that  changes  in  Grashof  number  and  in¬ 
clination  angle  have  a  stronger  influence  on  heat  transfer 
in  a  rectangular  enclosure  where  A=1.0  than  in  enclosures 
with  A=2.0  or  3.0.  This  is  due  to  the  side  wall  distur¬ 
bances  in  a  square  enclosure  occupying  a  larger  proportion 
of  the  enclosure  than  in  higher  aspect  ratio  cavities. 

In  summary,  the  numerical  results  for  the  inclined  mod¬ 
erate  aspect  ratio  rectangular  enclosure  indicate  the  val¬ 
idity  of  the  numerical  procedure  and  computer  code. 


Effect  of  Partitions 

In  this  section,  the  validity  of  natural  convection 
heat  transfer  results  for  composite  rectangular  enclosures 
will  be  determined.  Then,  the  heat  transfer  results  for 
the  partitioned  enclosure  will  be  compared  to  heat  trans¬ 
fer  results  of  rectangular  enclosures  without  partitions. 
From  these  comparisons  the  effectiveness  of  partitions  in 
reducing  natural  convection  heat  transfer  across  and  en¬ 
closure  will  be  resolved. 

Composite  Enclosures.  To  determine  the  validity  of 
composite  heat  transfer  solutions,  two  partitioned  rectangu¬ 
lar  enclosures  were  investigated.  These  enclosures  are 
shown  in  Figure  (9).  The  square  enclosure  is  composed  of 
two  triangular  cells  while  the  rectangular  enclosure  (A=2) 
is  composed  of  three  triangular  cells.  Natural  convec¬ 
tion  heat  transfer  coefficients  were  evaluated  for  each  of 
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these  triangular  cavities.  The  heat  transfer  coefficients 
of  adjoining  triangular  cells  were  then  compared  at  the 
common  partition  boundary. 

In  the  square  enclosures,  the  total  heat  transfer  a- 
cross  the  partition  was  balanced,  however,  heat  conduction 
along  the  partition  was  required  to  maintain  the  partition 
boundary  linear  temperature  distribution.  On  the  other 
hand,  the  total  heat  transfer  across  the  partitions  in  the 
rectangular  enclosure  (A=2.0)  did  not  balance,  eventhough 
the  total  heat  transfer  at  the  hot  and  cold  walls  were  bal¬ 
anced.  These  results  indicate  valid  composite  heat  trans¬ 
fer  solutions  exist  only  for  square  enclosures  with  a  diagon¬ 
al  partition.  An  iterative  numerical  procedure  is  required 
for  other  than  square  enclosure. 

Partitioned  Square  Enclosures.  To  determine  the  effec¬ 
tiveness  of  diagonal  partitions  in  suppressing  natural  con¬ 
vection,  heat  transfer  coefficients  were  evaluated  for  three 
square  enclosures  with  adiabatic  side  walls  at  various  in¬ 
clination  angles,  and  Grashof  numbers.  The  results  are 
shown  in  Figure  10. 

From  this  figure,  it  was  concluded  that  a  diagonal  par¬ 
tition  in  a  square  enclosure  reduces  natural  convection  heat 
transfer  by  more  than  50%.  An  explanation  for  the  substan¬ 
tial  decrease  in  natural  convection  heat  transfer  across  the 
square  enclosure  is  that  the  diagonal  partition  increases 


the  wetted  surface  area  by  70%.  And,  when  air  flows  over 
this  surface,  additional  friction  forces  are  generated  which 
suppress  the  natural  convection  in  the  enclosure. 

Figure  10  also  illustrates  the  effect  of  partition  angle 
on  natural  convection  heat  transfer.  The  partition  oriented 
at  $=-45  degrees  to  the  hot  wall  provided  an  additional  15% 
decrease  in  natural  convection  heat  transfer  when  compared  to 
the  partition  at  <f>=45  degrees  .  To  provide  an  expanation 
for  this  further  decrease  in  heat  transfer,  the  steady  state 
temperature  distributions  for  the  three  square  enclosures 
were  examined  In  Figures  11,  12  and  13. 

From  these  figures,  it  was  concluded  that  the  partition 
angle  affects  the  thickness  of  the  thermal  boundary  layer  at 
the  hot  wall.  The  partition  oriented  at  $=-45  degrees  re¬ 
sulted  in  the  greater  reduction  in  heat  transfer  because  it 
produced  the  thicker  thermal  boundary  layer  on  the  lower  half 
of  the  hot  wall  (the  region  of  highest  heat  transfer),  as 
evidenced  by  the  increased  spacing  between  isotherms. 

In  summary,  the  composite  heat  transfer  solution  was  val 
idated  for  square  enclosures  with  partitions.  Then,  a  diagon 
al  partition  in  an  inclined  enclosure  was  determined  to  re¬ 
duce  natural  convection  heat  transfer  by  more  than  50%. 
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V.  Conclusions 


Conclusions 

\ 

The  results  of  this  study  provide  the  following  con¬ 
clusions  : 

1.  The  finite  difference  numerical  procedure  yields 

valid  transient  and  steady  state,  natural  convection  heat 
transfer  coefficients  for  air  filled,  inclined,  moderate 
aspect  ratio,  rectangular  enclosures.  The  procedure  is 
limited  to  GR<2.0xl05  ,  and  lO^S^-90 

2.  A  thin  diagonal  partition  in  an  inclined  square  en¬ 
closure  with  insulated  side  walls  reduces  natural  convec¬ 
tion  across  the  enclosure  by  more  than  50%. 


VI.  Recommendations 


Recommendat ions 

The  following  recommendations  are  suggested  for  fol¬ 
low  on  investigations: 

1.  Numerical  techniques  should  be  developed  which 
would  extend  the  present  procedure's  Grashof  number  limi¬ 
tation  of  2.0x10s  . 

2.  A  numerical  procedure  should  be  developed  to  eval¬ 
uate  natural  convection  heat  transfer  in  an  inclined  rec¬ 
tangular  enclosure  with  more  than  one  partition.  One  ap¬ 
proach  to  this  procedure  would  be  to  assume  an  initial  par¬ 
tition  temperature  distribution,  compute  the  flow  field  in 
each  partition  cell,  then  compare  the  resulting  heat  trans¬ 
fer  coefficients  at  each  partition  grid  point.  If  the  heat 
transfer  coefficients  do  not  balance,  select  a  new  partition 
temperature  distribution,  and  repeat  the  process  until  the 
heat  transfer  coefficients  balance  at  each  partition  grid 
point.  This  procedure  preserves  the  computational  effi¬ 
ciency  obtained  with  tridiagonal  matrices. 

3.  The  effects  of  radiation  heat  transfer  within  the 
enclosure  on  partition,  and  side  wall  thermal  boundary  con¬ 
ditions  should  be  investigated. 
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Appendix  A 


Numerical  Data 


Triangular  Cell  with  <p 


DX  =  D 


0.1  Dt 


gr 

AR 

s 

Nu 

20000 

1.0 

90 

1.5088 

60 

1.5885 

1.0 

30 

1.4610 

60000 

1.0 

90 

60 

2.7108 

30 

2.5536 

100000 

1.0 

90 

60 

3.389 

30 

2.5536 

200000 

1.0 

90 
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4.5353 

30 

3.7664 

20000 

1.0 
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1.0 
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2.61767 

90 

2.9630 

60 

3 . 6636 

30 

3. 7764 

200000 


1.0 


Rectangular  Cavity 


JJlWJJ  i  Hi 


Trrrmmrm 

DX  =  DY  -  0.1  DT  =  .002 


. gr 

AR 
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Nu 
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3.036 
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30 

6.2863 
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Appendix  B 


I.  Inclined  Rectangular  Enclosure  with 

Insulated  Side- 

walls. 

II.  Inclined  Triangular  Enclosure  $=45 

degrees , 

In- 

sulated  Side  Wall. 

III.  Inclined  Triangular  Enclosure  $=45 

degrees , 

In- 

sulated  Side  Wall. 

IV.  Inclined  Triangular  Enclosure  Two  Partition 

Walls. 

V.  Tridiagonal  Algorithm  Subroutine. 

Computer  Code  Notation 

AR 

Aspect  Ratio 

DT 

Time  Increment 

DX 

Grid  Spacing  X-direction 

DY 

Grid  Spacing  Y-direction 

GR 

Grashof  Number 

IX 

Number  Grid  Points  in  X-direction 

IY 

Number  Grid  Points  in  Y-direction 

PR 

Prandil  Number 

S 

Inclination  Angle  (Degrees) 

SI 

Inclination  Angle  (Radians) 

ST 

Stream  Function 

T 

Dimensionless  Temperature 

U 

Dimensionless  Velocity  in  X-direction 
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Thomas  K.  Toltzien  was  born  on  May  15,  1949  in  Madison, 
Wisconsin  to  Paul  R.  and  Jo  K.  Toltzien.  In  June  1967,  he 

graduated  from  Madison  West  High  School.  He  then  attended 
the  University  of  Wisconsin-Madison ,  and  received  a  Bachelor 
of  Science  degree  in  Mechanical  Engineering  in  January  1972. 
He  entered  the  United  States  Air  Force  in  June  1972,  and  re¬ 
ceived  a  commission  through  OTS  in  September  1972.  He  was 
then  assigned  to  Mather  AFB,  Ca.  for  navigator  training. 
Following  graduation  in  July  1973,  he  was  assigned  to  the 
16th  Tactical  Reconnaissance  Squadron  (TRS),  Shaw  AFB,  SC 
for  RF-4C  Weapon  Systems  Officer  (WSO)  training.  Upon  com¬ 
pletion  of  training,  he  was  assigned  to  the  14th  TRS,  Udorn 
RTAFB,  Thailand  as  a  Mission  Ready  RF-4C  WSO.  In  April  1975, 
he  was  reassigned  to  the  1st  TRS,  RAF  Alconbury,  England. 
While  stationed  at  RAF  Alconbury,  he  was  a  RF-4C  Instructor 
Weapon  Systems  Officer,  and  a  Flight  Examiner.  Capt .  Tolt¬ 
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